UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

LIBRARY ROUTIN E M22 - 2^1 
By C W. Gear 

TITLE Eigenvalues and Eigenvectors of a Symmetric Matrix Using 

the Methods of Givens and Wilkinson. (SADOI Only) 
TYPE Complete Program 

NUMBER OF WORDS 10^4 + data storage 
MAXIMUM SIZE OF MATRIX 65 x 65 



TIME 



ACCURACY 



PROGRAMMED CHECKS 



Time for all roots and vectors is a maximum of 

.096n + .00 764n 2 + °0Q0 213n 5 mins. 

+ read in time. Where n is the order of the matrix. 

If only p values and q > p vectors are found the time is 

approximately 

p[.096 + .00237 n] + q .0.0014 n 

+ .00 287 n 2 + .000 213 r? mins. 
If the eigenvalues are distinct, then the vectors are 
good to ten places or better. The apparent accuracy 
will increase with n as the vectors are normalised [to 
length 1] , and hence as n increases the relative size of 
each decreases . If any eigenvalues are close or equal, then 
the eigenvectors will not be orthogonal for those close roots, 
but in all cases tried, appear to span the subspace generated 
by the actual eigenvectors. For further details see the 
description. 
FF 03F in location 209 (sexadecimal) Incorrect no. of 

elements in matrix. 
FF 03L 173 Drum transfer error 

FF 040 213 

or 21F 

or 22F 



FF 0^1 



I83 



Sum of squares of 
elements in matrix 
exceeds 1. 
Sequence of 0's, l's, 
and 2's specifying which 
roots are to be found 
is in error, 
(see input) 



- 2 <- 

FF 0^2 1SK Drum transfer error 

FF 0^3 300 input of master routine 

in error 
INPOT (Size of matrix must not exceed 65) 

The master routine M22 is input in the usual fashion, If 
no input error occurs it will stop on 34 012 (sexadecimal). 
If an input error occurs it will stop on FF 0^3. In the 
former case, the data punched as below should be input with 
the black switch. 

The elements of the matrix should appear first, punched 
in the following orders 

a ll a l2 a 13 ' ' ' • a ln N 

a 22 a 23 a 2n N 

a 33 ' ' a 3n N 



a N 
nn 

After each row of the matrix an "N" or any of the characters 
"J", ,f F H , or "L" must appear. The elements must be punched 
as a sign followed by up to 12 decimal digits. The sum of 
the squares of all of the elements of the matrix must not 
exceed 1, or the machine will stop. Following the matrix 
a sequence of O's, l f s, and 2's must be punched, n in 
number, where n is the order of the matrix. They have the 
following significances 

If the eigenvalues of the matrix are ordered algebraically, 
the largest being the first, then if the kth digit iss 

Neither the kth root nor its corresponding 
vector will be calculated. 

1 The kth eigenvalue only will be calculated. 

2 Both the kth eigenvector and value will be 
calculated . 



OUTPUT 



- 3 - 



Following this two sexadecimal characters must be punched, 
each between 2 and N, the first specifies the number of 
decimal digits punched for the eigenvalues , the second 
those of the vectors . 

It will compute the required results and finally stop 
on 3k 015 on the R.H.S. of location OSN. A new matrix 
may be read in with the black switch. 
Example ; A tape might be punched as follows t 

+2 .+ +k -2 N +2 -2 +k N +2 + N +2 N (Matrix) 

2 2 1 3 5 (Sequence and print digits) 

Jor the above example the format would be as follows? 

+000000000007 (This is the difference of the sum of 



+800 



squares of the elements of the matrix 
before and after rotation. If this is 
unduly large, there may have been an 
arithmetic error) 



+50000 

-50000 

+50000 

-50000 

3 +000 

-50000 
-50000 

+50000 
+50000 

k -i+oo 



METHOD 



The first eigenvalue and vector were asked for, as were 

the third. The fourth eigenvalue only was requested. 

The roots are numbered. The eigenvectors are normalized. 

The method for finding eigenvalues is the same as in M20, 

except that the transformation matrix in the reduction to 

Jacobi form is calculated as well. 

If A is the original matrix, and J the Jacobi form, Q is 

defined by Q' A Q = J. 

If y is an eigenvector of the matrix J then x = Qy is 

easily seen to be an eigenvector of A. 



* k 



It [X^] are the roots of |j - XI | = and if X* is 

the machine approximation to the root X, , then (j *■ X* l) 

n 

is nearly singular. If b is any vector , say b = Z b. y 

1 X i 

where the y. are the eigenvectors of the matrix J, then 

the solution of (j - VyOy = "b • « • • • • • • (1 ) 

is y » (J - X* k l) S) - Z x i 

i=l x» - x*» 

Providing that — — — — dominates all other coefficients, 



then y is a good approximation to y, . This will happen if b, 
is not small, unless there should be two or more close roots. 
In this case, the apparent eigenvector will be in the space 
spanned by the eigenvectors corresponding to the adjacent 

roots. In this method it has been found that with roots 

-39 
different only by 2 ^ appreciably different vectors are 

found. They may not be orthogonal, but they define the 

space within which the true vectors lie. If two roots are 

equal in the machine then they are made different by 2 

in order to get different vectors . 

The solutions of equation (l) are found by triangularizing 

the matrix (j - X,l) and then back substituting „ The vector 

b is such that it consists of all l*s after the triangu- 

larization process has been carried out. The solution is 

then normalized, and premultiplied by Q to get the 

eigenvector of A. 



Igr 



DAT E April 18, 1958RT:$/20/ 59 
PROGRAMMED BYj£u>. 
APPROVED BY 
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NOTES 



PAGE 5 



M 22 



LOCATION 


01DER 


NOTES PAGE 6 


203 


42 205F 
L5 ^55F 






204 


J2 3F 

50 204F 






205 


26 (Pl6) 
L5 12F 






206 


00 IF 
kO F 






207 


92 967F 
L5 F 






208 


50 F 
50 208F 






209 


26 (Pl6) 
92 139F 






210 


26 187F 
L5 35F 






211 


ko 453F 
4i 45OF 






212 


L5 89^F 
LO 12F 






213 


40 63^F 
L5 959F 


j 




2l4 


4o 699F 


i 




kl j6kF 






215 


F5 450F 
40 450F 






216 


l4 46f 
42 22 OF 




* 


217 


F4 399F 
42 257F 






218 


L4 428F 




i 




42 221F j 




219 


L5 450F 








LO 34F 







M 22 



LOCATION 



220 
221 
222 
223 
224 
225 
226 

22? 
228 
229 

230 
231 



233 
234 

235 
236 



ORDER 



NOTES 



PAGE 7 



M 22 



36 264f 
L5 F 

40 IF 
L5 F 
40 2F 
L7 IF 
12 2F 
32 244F 
L5 450F 
l4 46f 

42 240F 
L4 31F 
42 239F 
L4 44f 
42 238F 
42 24lF 
L4 31P 
42 236F 
L4 44F 
42 243F 
f4 399F 
42 232F 
L4 32F 
42 233F 
LI 12F 
L4 F 
40 3F 
50 F 
L5 233F 
42 242F 
71 IF 
66 2F 
S5 F 
40 F 



LOCATION 


ORDER 


NOTES PAGE 8 


257 


50 3F 
Tl IF 






258 


66 2F 
L5 F 






239 


S4 F 
40 F 






240 


L5 2F 

40 F 






24l 


L5 5F 

40 F 






242 


22 242F 
L5 F 






2^3 


22 2^3F 
40 F 






244 


22 257F 
40 5F 






245 


L5 450F 
F4 46f 






246 


42 255F 

l4 44f 






247 


42 252F 
42 258F 






248 


L4 31F 
42 256F 






249 


F4 428F 








42 255F 






250 


L4 J2F 
42 254F 






251 


L5 3F 
52 258F 






252 


50 2F 
71 F 






255 


66 IF 
L5 F 







M 22 



LOCATION 



ORDER 



NOTES 



PAGE 9 



M 22 



254 


S4 F 




50 F 


255 


LO 12F 




ho F 


256 


S5 F 




40 F 


257 


22 257F 




41 F 


258 


26 215F 




L5 F 


259 


40 3F 




L5 2F 


260 


L4 IF 




40 4F 


261 


L3 4F 




36 263F 


262 


LI 3F 




40 3F 


263 


50 3F 




22 253F 


26k 


L5 34F 




L4 46F 


265 


42 266F 




FO 399F 


266 


42 268F 




L5 F 


267 


40 457F 




26 268F 


268 


L5 27F 




40 F 


269 


L5 450F 




LO 31F 


270 


40 450F 




L4 349F 



LOCATION 


ORDER 


NOTES PAGE 10 


271 


42 330F 
iA 31F 






272 


k2 283F 
iA 399F 






273 


42 3l4F 
42 280F 






2jk 


iA 32F 
42 282F 






275 


IA 32F 
42 331F 






276 


50 F 
42 287F 






277 


F5 283F 
k2 297F 






278 


F5 331F 
42 299F 






279 


4l 452F 
L3 450F 






280 


36 337F 
L3 F 






281 


36 kkjF 
50 F 






282 


22 282F 
50 F 






283 


22 283F 
71 F 






284 


40 IF 

36 290F 






285 


L5 457F 
36 291F 






286 


iA IF 
36 293F 






287 


4o 456F 

L5 F 







M 22 



LOCATION 


1 

ORDER 


NOTES PAGE 11 


288 


kO IF 
L3 IF 






289 


32 3l^F 
22 297F 






290 


L5 ^57F 
36 292F 






291 


Lk IF 
26 287F 






292 


Lk IF 
36 287F 






293 


10 IF 
Ik klF 






29^ 


ko k%F 

F5 ^52F 






295 


kO k^ZF 
L5 ^57F 






296 


10 IF 
kO ^7F 






297 


22 287F 
50 F 






298 


71 IF 
40 IF 






299 


L5 ^52F 
iA F 






300 


42 302F 
kO 2F 






301 


L3 2F 
32 303F 






302 


L5 IF 
10 F 






303 


40 IF 
L5 IF 






30k 


36 307F 

L5 456F 

■ '■' — t 







M 22 



LOCATION 


ORDER 


NOTES PAGE 12 


305 


36 308F 
Lk IF 






306 


32 310F 
22 308F 






307 


L5 k%¥ 
32 309F 






308 


lA IF 

ko 456F 






309 


22 31^F 
lA IF 






310 


32 308F 
10 IF 






311 


L^ klF 

ko k%j 






312 


F5 ^52F 
40 *I-52F 






313 


F5 ^57F 
10 IF 






31^ 


ko ij-57F 
L5 F 






315 


40 IF 
L7 IF 






316 


L2 456F 

50 4op 






317 


32 327F 
iH 3F 






318 


L5 ^56f 
36 326F 






319 


LI IF 
ko IF 






320 


F5 3F 
k2 3F 






321 


L7 ^56F 








10 IF 







M 22 



LOCATION 


ORDER 


NOTES PAGE 13 


322 


kO h^GF 
L2 IF 






323 


LO 31F 
36 320F 






32^ 


L5 3F 
k2 326F 






325 


lA ^52F 
kO k52F 






326 


L5 ^57F 
10 IF 






327 


ko ^57F 
L7 ^56F 






328 


L2 IF 
32 332F 






329 


L5 ^F 
66 if 






330 


S5 F 
40 F 






331 


L5 ^52F 
40 F 






332 


26 269F 
L5 ^56f 






333 


lA IF 
^0 3F 






33^ 


L3 3F 

36 336F 






335 


L5 27F 
22 330F 






336 


LI 27F 
22 330F 






337 


L5 3^F 
kO lt-50F 






338 


4l 10F 
kl 11F 







M 22 



LOCATION 


ORDER 


NOTES PAGE lk 


559 


kl ^52F 
L5 ^50F 






5^0 


iA 5^9F 
k2 550F 






5^1 


p4 ij.i4.3F 
42 5^2F 






5^2 


50 40F 
L5 F 






5^5 


00 IF 
\2 5^8F 






3kk 


40 IF 
LI IF 






5^5 


iM- kk9¥ 
52 546F 






5^6 


27 5^9F 
L5 IF 






5^7 


52 350F 
50 11F 






5^8 


L5 10F 
10 F 






5^9 


kO 10F 
S5 502F 






550 


kO 1XF 
L5 F 






551 


kO IF 
50 IF 






552 


75 IF 
kO F 






555 


L5 452F 
k2 557F 






55^ 


LI 452F 

lh kk9F 






555 


56 556F 
27 558F 







M 22 



LOCATION 


OBDER 


NOTES PAGE 15 


356 


L3 ^52F 
32 358F 






357 


L5 F 








10 F 




,. 


358 


40 F 
L5 11F 






359 


Sk F 
36 362F 






360 


Lk klF 

^0 11F 






361 


F5 F 
40 F 






362 


L5 F 
L^ 10F 






363 


kO 10F 
36 368F 






36k 


50 11F 

10 IF 






365 


L4 hlF 
kO 10F 






366 


S5 F 
^0 IIP 






367 


F5 ^52F 
40 452F 






368 


L5 J +50F 
LO 31F 






369 


40 ^50F 
LO 31F 






370 


32 339F 
50 il-OF 






371 


L5 ^52F 
10 IF 






372 


^0 ^52F 
S5 F 







M 22 



LOCATION 


ORDER 


NOTES PAGE l6 


373 


LO 51F 
52 597F 






37k 


L5 12F 
kO kF 






375 


L5 11F 
ko F 






376 


L5 lOF 
50 576F 






377 


26 (Rl) 
kO lOF 






378 


F5 450F 
ko 450F 






379 


Lk 177F 
k2 581F 






580 


k2 592F 
Fk kk^F 






581 


k2 590F 
L5 F 






582 


kO F 
L3 ^52F 






585 


56 595F 
50 ^OF 






38k 


L5 ^52F 
k2 586F 




• 


585 


LO 48F 
52 590F 






586 


L5 F 
10 F 






587 


kO F 
L7 F 






588 


L2 10F 
52 396F 






389 


L5 F 








50 il-OF 







M 22 



LOCATION 


ORDER 


NOTES PAGE 17 ; 


390 


66 10F 








L5 F 






391 


iA ^52F 
k-0 k^2F 




i 


392 


S5 F 
40 F 






393 


Fl ^50F 
iA 3^F 






39^ 


36 378F 
26 ^02F 






395 


50 *l-OF 
L5 F 






396 


26 387F 
50 27F 






397 


22 390F 
L5 10F 






398 


50 11F 
10 IF 






399 


40 10F 








S5 130F 






^00 


kO 11F 
F5 ^52F 






j 401 


40 J+52F 
26 37^ 






1 402 

j 

1 


L5 31F 
40 45^F 






403 


L5 ^F 
00 IF 






1 kok 


40 kF 
L5 ^55F 






405 


J2 3F 
50 405F 






ho6 


26 (Pl6) 
92 967F 







M 22 



LOCATION 



ORDER 



NOTES 



PAGE 18 



M 22 



407 


L5 4f 




22 4o8f 


4o8 


50 F 




50 4o8f 


409 


26 (pi6 




92 135F 


410 


L5 503F 




22 435F 


4ii 


F5 4^4f 




4o 454f 


412 


l4 444f 




4o 4l5F 


4l3 


F5 46f 




42 4i6f 


4l4 


42 420F 




26 415F 


415 


00 F 




00 F 


4l6 


22 4l6F 




40 F 


417 


L5 4l5F 




L4 32F 


4l8 


40 415F 




F5 4i6f 


419 


42 4i6f 




LO 445F 


420 


36 4i5F 




Ll F 


421 


L4 453F 




4o 453F 


422 


F5 42 OF 




42 420F 


423 


LO 446f 




32 420F 



LOCATION 


ORDER 


NOTES PAGE 19 


k2k 


L5 31F 
40 450F 






k25 


4l 10F 
50 40F 






k26 


L5 ^50F 
L4 46f 






k27 


42 429F 
LO 399F 






428 


42 430F 
S5 19^F 






429 


22 429F 
50 F 






430 


22 430F 
?4 F 






431 


Lk lOF 
40 lOF 






432 


F5 450F 
40 450F 






433 


LO 34F 
32 434f 






434 


26 426F 
L5 lOF 






435 


50 F 
50 435F 






436 


26 (Pl6) 
92 131F 






437 


92 515.F 
L5 fc^F 






438 


LO 34F 
32 439F 






1^39 


26 4llF 
L5 455F 






44o 


40 451F 
92 135F 







M 22 



LOCATION 


ORDER 


NOTES PAGE 20 


kkl 


L3 ^53F 
36 H3F 






kk2 


FF 66f 
FF 66f 






kk3 


26 I87F 
00 260F 






kkk 


85 11F 

00 6783F 






Kk<j 


K2 ^i6f 
ifO F 






kk6 


S6 415F 
LI F 






kkj 


L5 ^8F 
ho 4p2F 






kkQ 


L5 27F 
22 330F 






kk9 


00 F 
00 60F 






450 


00 F 
00 F 






k^l 


00 F 
00 F 






^52 


00 F 
00 F 






453 


00 F 
00 F 






.k-5k 


00 F 
00 F 






^55 


00 F 
00 F 






456 


00 F 
00 F 






457 


00 F 
00 F 







M 22 



LOCATION 


ORDER 


NOTES PAGE 21 




26 15N 




Interlude to store B on drum. 




00 24k 






24 


50 154F 








50 L 




Read routine B from drum 


25 


26 (Yl) 
00 IIOIOF 






26 


00 300F 

26 154f 
00 154k 






154 


4i 451F 




J = 




F5 451F 




J -> J + 1 


155 


4o 451F 








FO 3^F 




J = N + 1 ? 


156 


36 24F 
L5 451F 






157 


L4 197F 








42 158F 




Are we to find this (j) value? 


158 


32 158F 








Ll F 




[Ap - calculate eigenvalues from 


159 


32 154f 
L5 451F 




reduced form] 


160 


l4 46f 
42 163F 






161 


42 167F 
L4 32F 






162 


42 165F 
42 169F 






163 


L5 313F 
L4 F 






164 


10 IF 
40 F 






165 


50 296F 
Ll F 







M 22 



LOCATION 


ORDER 


NOTES 


PAGE 22 


166 
167 


10 IF 
Lk F 
32 15^F 
L5 F 




-38 
L T - U T + 2 ^ > ? 




168 


10 IF 
kO F 








I69 


50 297F 
F5 F 




u = (Uj + Lj)/2. 




170 


10 IF 
Lk F 








171 


ko 38F 
kl 450F 




1 = 




172 


kl ^52F 








173 


kl 3F 
in 8F 
L5 27F 




m , = 




17^ 


kO kF 
F5 450F 




m = 1 - 2~ 59 




175 


k2 450F 
FO 3^F 




I -*I + 1 




176 

1 


32 240F 
L5 ^50F 




I = N + 1 ? 




j 177 


Lk 192F 
42 188F 








178 


L4 32F 
k2 195F 




Set I dependent addresses . 




179 


42 212F 
42 233F 








180 


L4 32F 
42 183F 








181 


LI 38F 
40 7F 








182 


kl 6f 
L5 255F 




Link for scaling routine. 
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185 


hO 298F 
L5 F 




Scale a. 


184 


22 26kF 








L3 3F 




m i-2 " ° ? 


185 


32 223F 








L3 4f 




"i-1 " ° * 


186 


32 233F 








L5 256F 




Link for scaling routine. 


187 


kO 298F 








L5 8F 




A. = V. , + 6. V 4 t 
1 i-1 i-Zi-1 


188 


Lh 2F 
LO F 






189 


40 IF 








32 210F 




A. > ? 

1 — 


190 


L5 254F 








k2 200F 




It, -» Jf . 
1 J 


191 


50 F 








75 4f 




P i ^-1 


192 


40 7F 
S5 763F 






193 


kO IIF 
LI IF 






19^ 


^2 196F 
LO kSF 






195 


32 209F 








50 F 




A > 40 ? 


196 


75 3F 








10 F 




q i-l ' m i-2 


197 


kO 5F 








S5 568F 




Subtract 


198 


ho 6f 
L5 2F 






199 

\ 


40 8F 
L5 IIF 
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200 


26 257F 
26 F 




Scale to m and E. 




201 


L5 2F 
IX 8F 









202 


ko 8f 

L5 ^F 




v. = E. + 0, 

1 i i 




205 


kO 5F 
L5 F 








26k 


kO kF 
L5 5F 




m. . > ? 
i-l — 




205 


56 208F 
L5 ^F 








206 


52 17^F 
F5 i i-52F 




m. > ? 

■ 1 — 

p -»p + 1 




207 


k2 452F 
22 17^F 








208 


L5 ^F 
52 206F 




m. > ? 
1 — 




209 


22 17^F 
kl 5F 








210 


26 I98F 
k2 219F 




*2 _M, j 




211 


L5 221F 
k2 200F 




4 i-X ' m i-2 




212 


50 5F 
75 F 




P i m i-1 




215 


kO 5F 
S5 F 








21^ 


ko 6f 

50 F 




-■ 




215 


75 ^F 
kO 7F 








216 


L5 IF 
52 220F 




A = ? 
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217 


Lk 48F 
36 219F 




A > kO ? 




218 


22 222F 
00 F 








219 


L5 7F 
10 F 








220 


ko 7F 
L5 2F 








221 


kO 8F 
S3 231F 








222 


26 257F 
kl 7F 




Scale to m. and Z. 
i 1 




223 


26 257F 

Li Af 




m > ? 




22*1- 


36 228F 
L5 27F 








225 


kO 3F 
L5 F 




m. . = 1 - 2~ 39 
l-l 




226 


ij-0 kF 
L5 2F 




m. = p. 
1 1 




227 


4o 8f 

22 204F 








228 


LI 27F 
kO 3F 




m. t = -1 + 2~ 59 
i-l 




229 


LI F 
kO kF 




m.~ = - P± 




230 


22 204F 
00 F 








231 


LI IF 
Lk 8F 








232 


L^ 2F 
ko 8F 








233 


22 202F 
L3 F 




*i-l " ° * 
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234 


32 224F 








L5 3F 




m i-2 > ° * 


235 


32 237F 








L5 27F 






236 


40 4f 








4l 3F 




in. = 1 - 2" 59 
1 


237 


22 204F 








LI 27F 




m. = -1 + 2" 59 

1 


238 


26 236F 

L5 4f 






239 


40 3F 

41 kF 






240 


22 204F 








L5 451F 


* 


xr _ t 
JX = O 


24l 


kO k^kF 








L5 454F 




' 


242 


L4 46f 








k2 251F 




Set K dependent 


243 


42 253F 
L4 52F 




addresses . 


244 


42 248F 
42 246f 






2k5 


L5 ^52F 
LO 454f 






246 


32 25IF 








Ll F 




~1U + u > ? 


2^7 


L4 38F 
36 249F 






248 


L5 38F 






249 


40 F 
F5 4^4f 




U K = u 




42 4^4f 




K ->K + 1 


250 


FO 34F 






i 


36 163F 




K = n + 1 ? 
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251 


22 2*kLF 
L5 F 




1^ - u > ? 




252 


LO 38F 
36 2^9F 








253 


L5 38F 
^•0 F 




L K = u. 




25b 


26 2^9F 
00 201F 








255 


22 18^F 
22 l8^F 




Constants for exit 




256 


22 200F 
22 238F 




from scale routine. 




257 


LO 6f 
32 262F 




Form least significant 




258 


L^ 4lF 










ho 6f 




half and test for 




259 


L7 5F 
36 26^F 




carry. This part 




260 


L^ 7F 
36 299F 




does a double precision 




26l 


LO 31F 
36 299F 




subtraction before 




262 


22 267F 
kO 6F 




scaling of 




263 


LI 5F 
22 26ii-F 




N(7, A) -N(5, 6) 




26^ 


Fl 5F 
ho 5F 




Scale routine 




265 


36 302F 
L5 7F 








266 


36 303F 
Lk 5F 




Number is made negative 




267 


36 299F 
kO F 


1 
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268 


4l 2F 
L5 169F 








269 


k6 296F 
L3 F 




Test for 0. 




270 


32 298F 

50 6f 




. 




271 


L5 F 
00 IF 


f-1 








272 


36 282F 












00 UF 






Scale by 2 




273 


36 283F 
00 IF 






until positive 




2jh 


32 2QkF 
00 IF 










275 


36 286F 
00 IF 










276 


32 287F 
00 IF 










277 


36 289F 
00 IF 










278 


32 290F 
00 IF 










279 


36 292F 
kO F 










280 


F5 265F 
Lk 2F 






Increase scale factor by 8 




281 


k2 2F 
26 27XF 










282 


k-0 F 
27 293F 








283 


kO F 
L5 31F 




Increase scale factor by 




284 


26 293F 
kO F 




appropriate amount 




285 


F5 31F 
26 293P 
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286 


kO F 
L5 236F 






287 


26 293F 
kO F 






288 


F5 236F 
26 293F 






289 


kO F 
L5 26^F 






290 


26 293F 
kO F 






291 


F5 264F 








26 293F 




; 


292 


kO F 
L5 265F 






293 


Lk 2F 








k2 2F 




Add to scale factor 


29^ 


L5 F 








10 IF 




Adjust scaled number to correct 


295 


iA klF 
kO F 




sign and -f by 2. 


296 


22 F 
LI F 






297 


ho F 








26 298F 




■ 


298 


22 F 








22 F 




exit. 


299 


10 IF 








Lk kUP 




Overflow adjustment, 


300 


kO F 
LI 31F 






301 


kO 2F 
26 298F 






302 


L5 7F 
32 30^F 
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303 


Lk 5F 






. 


36 306F 






30k 


22 267F 
L^ 5F 






305 


36 306f 
26 299F 






306 


kO F 
il-1 2F 






307 


L5 165F 
k6 296F 






308 


LI 6f 
32 3HF 






309 


Li*- 41F 
ij-0 6f 






310 


Fl F 








40 F 






311 


22 269F 
LI F 






312 


kO F 








22 269F 






313 


00 F 
00 2F 






31k 


41 35F 








26 508F 








00 K 






315 


L5 29F 








kO kL 




[A, s reduction of matrix to 


316 


Lk 32F 




Jacobi form and read in of 




kO 10L 




matrix] . 


317 


F5 i*-3F 
i+2 % 






318 


Fk 32F 








42 11L 




Read out from drum of reduced matrix* 
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319 


85 11F 








00 2560F 






320 


22 5L 








40 F 


: 




321 


L5 4l 








F4 32F 






322 


40 4L 








F5 5L 






323 


42 5L 


1 






LO 15L 






324 


36 i6l 

26 10L 






325 


85 11F 








00 2560F 

■ 






326 


22 11L 

40 F 


t 




327 


L5 10L 
F4 32F 






328 


40 10L 








F5 11L 


! 




329 


42 11L 

26 4l 






330 


22 5L 
4o 959F 






331 


4l 959F 
4l 392L 






332 


4l 10F 
kl 11F 






333 


F5 392L 
42 392L 






33^ 


FO 34F 
32 51L 






335 


L5 392L 






! 


L4 42F 




. 
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336 


FO 31F 








42 47L 






337 


L4 32F 








42 46L 






338 


L4 32F 




s 


.. 


42 27L 


: 


5 


339 


42 36L 






; '• 


42 37L 




; 


340 


L4 32F 








42 28L 






34l 


42 29L 
42 30L 






342 


41 3F 








LI F 






343 


L4 33F 








LO F 






344 


40 33F 








50 F 






345 


22 30L 








75 F 




Sum squares 


346 


00 IF 








40 F 






; 347 


S5 F 








40 IF 






348 


L4 11F 
10 39F 






3^9 


40 2F 








L7 2F 




350 


L4 10F 
40 10F 






351 


S5 F 
50 F 






352 


22 37L 
74 F 
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353 


Lk lOF 
lA F 






35^ 


kO lOF 








S5 959F 






355 


kO 11F 
L3 F 


i 


i 


356 


32 k8h 


1 

: 






50 IF 




: 


357 


L5 F 








00 IF 






358 


kO F 








F5 3F 


^ 


I 


359 


42 3F 
L5 F 






360 


32 42L 








10 IF 






361 


Lk 4if 








40 F 






362 


L5 3F 
40 F 






363 


26 18L 
L3 IF 






364 


36 50L 








22 4lL 






365 


L5 28F 








k-0 3F 


: 




366 


23 46l 
92 139F 






367 


L5 37F 
LO 10F 






368 


50 12F 








50 53L 




Print difference of sum of 


369 


26 (Pl6) 
92 135F 




squares 
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370 


L3 33F 








36 57L 




Sum check on drum for A. 


371 


FF 63F 








00 568F 






372 


kl 39F 








F5 39F 






373 


kO 39F 








FO 3kF 






37^ 


36 65L 








L5 39F 




' . 


375 


Lk 6kL 




; ; 




k2 62L 




. 


376 


IX 32F 
k2 63L 




■ ' 


377 


LI 38F 
kO 6^kF 






378 


L5 38F 
14-0 699F 






379 


22 57L 

00 633F 






380 


in 39F 
F5 39F 






381 


k2 39F 
FO 3kF 






382 


32 73L 
L5 39F 






383 


IX 56L 








42 69L 




Read in sequence of 0*s, l's, 


38fc . 


81 k-F 




and 2 ' s . 




kO F 






385 


kO F 
LI F 






386 


Fk 31F 
36 73L 
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387 


FF 65F 








FF 65F 






388 


22 65L 




■ 




81 4f 






389 


00 20F 








46 9F 




Read in digits specifying no. 


390 


81 4f 

00 20F 




of places to print. 


391 


46 10F 






' 


26 154f 






392 


40 F 




Drum subroutine 




K5 102L 




II 


393 


42 82L 
L4 31? 






394 


42 101L 
10 20F 






395 


42 94L 
42 98L 






396 


42 102L 
42 107L 






397 


22 82L 

l4 f 






398 


42 111L 
42 112L 






399 


42 113L 
42 114L 






400 


S5 F 
36 88L 






401 


L5 77L 
42 92L 






402 


L5 116L 
22 89L 






403 


L5 115L 
42 92L 
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4o4 


L5 117L 








L4 F 






405 


50 40F 
00 6f 






4o6 


l4 f 
40 93L 






407 


40 103L 
26 F 






4o8 


00 F 
00 F 






409 


32 94L 

40 F 






4l0 


F5 93L 
40 93L 






4n 


F5 94l 
42 94l 






4l2 


LO 111L 
36 93L 






4l3 


L5 35F 
LO F 






4l4 


40 35F 
F5 98l 






415 


42 98L 
LO 112L 






4i6 


36 98L 
26 F 




■ 


4l7 


32 102L 

L5 F 






4l8 


00 F 
00 F 






4l9 


F5 103L 
40 103L 






420 


F5 102L 
42 102L 
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431 


LO 113L 
32 102L 






422 


L5 35F 
L4 F 






423 


40 35F 
F5 107L 






424 


42 107L 
LO ii4l 






425 


36 107L 
22 101L 






426 


S2 94l 

40 F 






427 


75 35F 

LO F 






428 


S2 102L 
L5 F 






429 


75 35F 








L4 f 




Constants 


430 


00 F 
00 93L 






431 


02 384F 
2N 105F 






432 


02 320F 
2N 105F 






433 


40 3F 




Drum subroutine 




K5 F 




I 


434 


42 126L 
10 20F 






435 


42 2F 
36 123L 






436 


L5 182L 
42 148L 






437 


F5 181L 
22 124L 
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438 


L5 140L 
42 l48L 






439 


L5 181L 
42 137L 






Ho 


42 l43L 
iH IF 






44i 


43 139L 
L5 F 






kk2 


42 IF 
F5 126l 






443 


42 i66l 
L5 2F 






444 


L4 IF 
42 153L 






445 


42 I67L 
42 162L 






446 


42 17TL 
F5 IF 






kkj 


LO 3F 
36 l49L 






448 


L5 3F 
LO 31F 






449 


4o F 
L4 2F 






450 


42 I83L 
42 184L 






^51 


50 F 
75 32F 






I4-52 


S5 F 
L4 F 






453 


l4 IF 
4o 152L 






454 


4o 168L 
50 F 
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455 


L5 2F 
S4 152L 






456 


42 158L 
42 172L 






457 


75 32F 
S5 5F 






458 


L4 3F 
l4 f 






459 


LO 31P 
40 157L 






46o 


40 173L 
L5 2F 






46l 


L4 34f 
42 185L 






462 


42 186L 
42 187L 






463 


42 188L 
26 F 






464 


L5 31F 
42 139L 






465 


L5 148L 
L4 142L 






466 


42 148L 
22 139L 






467 


00 F 

00 F . 






468 


32 153L 
40 F 






469 


F5 152L 








40 152L 






470 


F5 153L 
42 153L 






471 


LO 183L 
36 152L 
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^72 


00 F 
00 F 






473 


52 158L 
40 F 






474 


L5 157L 
L4 52F 






475 


40 157L 
F5 158L 






A76 


42 158L 
LO I85L 






477 


56 157L 
LI F 






478 


L4 33F 
40 33F 






479 


F5 162L 
42 162L 






48o 


LO I87L 
32 162L 






481 


32 166L 
26 F 






-482 


32 I67L 
L5 F 






483 


00 F 
00 F 






484 


F5 168L 
4o 168L 






485 


F5 167L 
42 167L 






486 


LO 184L 

32 167L 






487 


32 172L 

L5 F 






488 


00 F 
00 F 
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489 


L5 173L 
L4 32F 






490 


40 173L 
F5 172L 






491 


42 172L 
LO 186L 






492 


32 172L 
L5 F 






493 


L4 33F 
40 33F 






k9k 


F5 177L 
42 177L 






495 


LO 188L 
32 177L 






496 


22 166L 
00 29F 






497 


00 F 

00 167L 






498 


S2 153L 
40 F 






499 


S2 167L 
L5 F 






500 


S2 I58L 
40 F 






501 


S2 172L 
L5 F 






502 


S6 157L 
LI F 






503 


S2 172L 
L5 F 






504 


L5 198L 
L4 36F 






505 


46 198L 








10 20F 
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506 


k-2 212L 
k2 21JL 








507 


k2 214L 
26 198L 








508 


L5 208L 
k6 198L 




Read in begins here. 




509 


L5 39L 
k2 212L 




Reset routine 




510 


k-2 213L 
1*2 214L 








511 


4l 10F 
kl 11F 








512 


43 209L 
kl 2klL 








515 


50 959F 
50 198L 








5l4 


26 (N12 ) 
iH F 




Read in 1 row- 




515 


L5 21(N12) 
26 409L 




Find no. of elements read in. 




516 


k2 3^-F 










L5 3^F 




On first row set loco 3^ - to 


N, 


517 


iA ^08L 

ko 238L 




on each subsequent row check 
to see that N-l> N-2, ... 1 




518 


L5 236L 
kO 201L 




elements are read* 




519 


26 207L 
LO 233L 








520 


1(0 F 

L3 F 






1 


521 


56 207L 
FF 62F 








522 


F5 21+1L 
k2 2kll> 
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523 


JO 959F 
50 208L 




Store row on drum 




524 


26 118L 
00 F 




using routine I 




525 


L5 24lL 
42 209L 








526 


4l 12F 
4l 13F 








527 


L5 39F 
L4 959F 








528 


4o 39F 
50 959F 




Sum squares. 




529 


L5 13F 
74 959F 








530 


10 IF 










32 2l6L 




Test to see if sum 


exceeds 1. 


531 


FF 64f 
40 12F 








532 


S5 959F 
40 13F 








533 


L5 239L 
40 215L 








534 


F5 213L 
42 213L 








535 


42 2l4L 
42 212L 








536 


LO 238L 
36 212L 








537 


L5 11F 
L4 13F 








538 


10 39F 
40 F 








539 


S5 F 










40 11F 1 
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5^0 


LI F 
L^ 10F 






5^1 


Lk 12F 








32 227L 






5^2 


FF 6kf 




. 




kO 10F 






5^3 


L5 240L 


> 


' 




kO 215L 






5kk 


L5 23^L 




., 




LO 31F 




: 


5^5 


hO 233L 




■ ; 




LO 235L 




Test for last row 


5^6 


36 189L 








L5 237L 






5^7 


ho 201L 

26 2^2L 






5^8 


00 F 








00 F 






5^9 


00 F 
00 F 






550 


LL ^095F 
LO IF 






551 


22 20^L 
L5 3^F 






552 


k2 3J4-F 








L5 3^F 




Constants 


553 


NO 39F 
50 F 






55^ 


Ik 12F 
32 216L 






555 


10 IF 
32 216L 






556 


00 F 
00 F 
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557 


L5 10F 
00 IF 






558 


36 24^L 




i 




FF 6kF 




; 


559 


kO 37F 




1 




50 2¥<-L 




Square root of sum of squares 


560 


26 (Rl) 

ko 38F 






561 


26 2^7L 
00 F 






562 


L5 257L 
kO 2k&L 






563 


87 11F 








00 6789F 




Set Q = I on drum 


%k 


F5 2^8L 
kO 2^8L 






565 


LO 259L 
36 2k&L 






566 


L5 258L 
kO 253L 






567 


22 252L 
L5 27F 






568 


86 11F 
00 6789F 






569 


L5 253L 
F^ 32F 






570 


kO 253L 
LO 260L 






571 


32 252L 
26 26lL 






572 


87 11F 

00 6785F 






573 


86 11F 

00 6785F 
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576 
577 
578 
579 
580 
581 
582 
585 
584 

585 
586 

587 
588 

589 
590 



ORDER 

07 11F 
00 10946F 
06 lip 
00 109^6F 

41 392L 
L5 34F 
LO 31F 

42 389L 
42 282L 
L5 392L 
42 279L 
42 291L 
42 377L 
42 386L 
P5 392L 
40 392L 
40 392L 
F5 392L 
LO 34F 
36 L 

L5 392L 
L4 43P 
42 296L 
L4 31F 
42 362L 
42 371L 
42 324L 
L4 44f 
42 297L 
42 374L 
L4 31F 
42 327L 
42 366L 
42 370L 
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i = 

[Begin rotation] 



Set test constants 



i -+i +1 



Test for end of rotation 
i = n-1 ? 



Set i dependent addresses. 



LOCATION 


ORDER 
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591 


k2 292L 
k2 378L 








592 


F5 392L 
kO 393L 




j = i+1 




593 


50 89^F 
50 278F 




Read out one row of A» 




59^ 


26 ll8L 
00 F 








595 


L5 392L 
22 281L 








596 


50 766F 
50 281L 




Read out one row of Q<, 




597 


26 77L 
00 F 








598 


F5 393L 
\2 391L 








599 


L5 393L 
k2 295L 




Set j -dependent addresses 




600 


k-2 381L 








601 


k2 36^L 
\2 369L 








602 


h2 326L 
Yk kkF 








603 


k-2 368L 
42 373L 








60k 


42 325L 
F5 393L 








605 


50 959F 
50 2.90L 








606 


26 118L 
00 F 




Read out a row of A 


i 
1 

> 

1 


607 


L5 33F 








_ . 


L4 F 




1 


i 
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608 


40 33F 
L5 393L 






609 


50 83OF 








50 29^L 




Read out a row of Q 


610 


26 77L 








00 F 






611 


22 296L 








L5 F 




: 


612 


ko 3F 








L5 F 






613 


kO k-F 








L3 4F 




s 


6lk 


32 375L 


* 
\ 


[ 

i: 
* 




L7 3F 


1 
1 




615 


L2 4f 


\ 


r 




36 398L 






6l6 


50 3F 




' 




L5 4f 




L 


617 


40 3F 




1 




S5 F 




1 


618 


40 4f 


■ 




: 


L5 397L 


i 


J 


619 


42 317L 








46 324L 






620 


50 40F 








L5 4f 






621 


66 3F 






622 


S5 F 
40 5F 
7J 5F 




Calculate C. . and S. . 


623 


10 IF 
40 IF 






624 


L3 IF 
56 407L 
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625 


LJ IF 








40 IF 


1 




626 


50 40F 








6l IF 


; 




627 


41 F 
K5 F 


\ 


! 


628 


36 314L 








26 407L 






629 


S5 F 




' 




50 314L 






630 


26 (Rl) 








40 IF 


■ 




631 


50 IF 








7J 5F 






632 


40 F 








40 F 






633 


50 F 
79 F 




, 


634 


40 F 
50 IF 






635 


79 IF 
L4 f 






636 


LO 4lF 








50 40F 






637 


66 IF 

S5 F 






638 


10 IF 
L4 IF 






639 


40 F 
L5 F 






64o 


22 325L 
L4 F 






64l 


40 2F 








L5 F 




in 4 
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642 


22 527L 






■ 


40 F 






643 


22 328L 








L5 392L 




Prepare for (i,j) rotation of A 


644 


40 394L 


; 


: 




L5 394L 






645 


L4 43F 








42 334L 


; 




646 


42 3^L 






. 


42 340L 


i 




647 


F4 44f 


! 


i 




42 33^L 


' 


'> 


648 


42 338L 








42 34XL 






649 


50 5F 


i 


; 




75 F 






650 


40 F 




i 




S5 F 




Kotate in (i„j)th plane 


651 


50 6F 








74 F 




' 


652 


L4 f 








40 IF 






653 


50 5F 
75 F 






654 


40 F 
S5 F 






655 


50 6f 

70 F 






656 


L4 F 
4o F 






657 


L5 IF 
40 F 






658 


F5 394L 




Test for finish of (i,j)th 




42 39^L 




rotation 
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659 


FO 3kF 
32 3^5L 






660 


22 529L 

L5 31F 






661 


^0 395L 








L5 ^2F 






662 


IA 395L 








h2 351L 




Prepare for (i,j) rotation of Q 


663 


k2 357L 








^2 359L 






66k 


iA hhlp 
h2 353L 






665 


42 355L 
k2 358L 






666, 


50 5F 
75 F 






667 


lj-0 F 
S5 F 






668 


50 6f 








jh F 






669 


Lk F 
kO IF 






670 


50 5F 
75 F 






671 


ho F 
S5 F 






672 


50 6f 








70 F 




Rotate . 


673 


iA f 
40 F 






67^ 


L5 IF 
li-0 F 






675 


F5 395L 
40 395L 
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676 


LO 391L 








32 3^6L 




Test for finish. 


677 


50 5F 








75 F 






678 


kO F 








S5 F 




Special case for (ii), (i j ) and 


679 


50 6f 

lk F 




(jj) elements of A« 


680 


Lk F 
kO IF 






681 


50 5F 
75 F 






682 


kO F 
S5 F 






683 


50 6f 








7k F 






6Qk 


Lk F 
40 F 






685 


22 370L 
kO F 






686 


L5 IF 
kO F 






687 


LI IF 
Lk 2F 






688 


22 373L 
kO F 






689 


22 37^L 
kl F 






690 


22 375L 
F5 393L 






691 


JO 959F 








50 376L 




Write jth row of A "back on drum. 


692 


26 118L 
00 F 
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693 


L5 33F 








LO F 


: 




69k 


40 33F 
L5 393L 


: 


• 


695 


JO 83OF 
50 38OL 




Write jth row of Q back on drum. 


696 


26 77L 
00 F 






697 


F5 393L 
ho 393L 






698 


LO 3^F 




z 




32 384L 




Test for j = n - 1 ? 


699 


26 283L 
F5 392L 






700 


JO 89^F 




'.; 




50 385L 




Write ith row of A back on drum. 


701 


26 118L 
00 F 






702 


22 387L 
L5 392L 






703 


JO 766F 








50 388L 




Write ith row of Q back on drum. 


704 


26 77L 




; 




00 F 






705 


22 263L 
00 F 






706 


80 F 








00 F 




Constants 


707 


00 F 
00 F 






708 


00 F 
00 F 






709 


00 F 






. 


00 F 
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710 


00 F 
00 F 






711 


00 5F 








00 6f 




Constants 


712 


00 6f 








00 5F 






715 


kO F 
L3 F 






Jlk 


32 400L 
L5 396L 






715 


26 30^L 
L5 k$F 






716 


kO 5F 








L5 ^F 




■ 


717 


LO 3F 
ll-O.F 






718 


L3 F 
32 405L 






719 


LI 45F 

ifO 6f 






720 


22 32^L 
L5 ^5F 






721 


40 6F 
22 32^L 






722 


L5 ^1F 
22 315L 






723 


NO 39F 
50 959F 






724 


LO 198L 
10 20F 






725 

1 


k2 F 
L5 F 






726 


Lk 235L 
LO 31F 
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727 



^0 23^L 
26 201L 

(N12 )00K 

OOK 
L3 F 
3k 18F 
FF 67F 
2k 18F 
L7 2148F 
74 269^F 
26 L 
26 IN 



Input routine. 



Sum check routine 



